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Abstract. We present a criterion that provides an easy sufficient condition in order that a coilection 
of Abeiian integrals has the Chebyshev property. This condition invoives the functions in the integrand 
of the Abehan integrals and can be checked, in many cases, in a purely algebraic way. By using this 
criterion, several known results are obtained in a shorter way and some new results, which could not be 
tackled by the known standard methods, can also be deduced. 



1 Introduction and statement of the result 

The second part of Hubert's 16th problem |15) asks about the maximum number and location of limit cycles of a 
planar polynomial vector fields of degree d. Solving this problem, even in the case d = 2, seems to be out of reach at 
the present state of knowledge (see the works of Ilyashenko 'VP and Li Jibin ^0] for a survey of the recent results on 
the subject). Our paper is concerned with a weaker version of this problem, the so-called infinitesimal Hubert's 16th 
problem, proposed by Arnold [T]. Let w be a real 1-form with polynomial coefficients of degree at most d. Consider 
a real polynomial H of degree d + 1 in the plane. A closed connected component of a level curve H = h is denoted 
by 7h and called an oval of H. These ovals form continuous families (see Figure [2]) and the infinitesimal Hilbert's 
16th problem is to find an upper bound V{d) of the number of real zeros of the Abelian integral 



The bound should be uniform with respect to the choice of the polynomial H, the family of ovals {7/1} and the 
form uj. It should depend on the degree d only. (In the literature an Abelian integral is usually the integral of a 
rational 1-form over a continuous family of algebraic ovals. Throughout the paper, by an abuse of language, we use 
the name Abelian integral also in case the functions are analytic.) 

Zeros of Abelian integrals are related to limit cycles in the following way. Consider a small deformation of a 
Hamiltonian vector field = Xh + eY, where 



Then, see |17l I2U] for details, the first approximation in e of the displacement function of the Poincare map of is 
given by ([T]) with uj — Pdy — Qdx. Hence the number of isolated zeros of 1(h), counted with multiplicities, provides 
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Xh = --Hyd^ + H^d^ and Y ^ Pd^ + Qdy. 
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an upper bound for the number of ovals of H that generate hmit cycles of for e « 0. The coefficients of P and Q 
are considered as parameters of the problem and so the function 1(h) splits as a linear combination 

oiolo{h) + aih{h) + . . . + q„_i/„_i(/i), 

where depends on the initial parameters and Ik{h) is an Abelian integral with either ui — x^y-'dx or uj — x^y-'dy. 
(In fact it is easy to see, using integration by parts, that only one type of these 1-forms needs to be considered.) 
Therefore the problem is equivalent to find an upper bound for the number of isolated zeros of any function belonging 
to the vector space generated by Ik{h) for fc = 0, 1, . . . , n — 1. This problem is strongly related to showing that the 
basis of the previous vector space is a Chebyshev system. In fact, the great majority of papers studying concrete 
problems on the subject show this kind of property. 

In this paper we focus on the case in which H has separated variables, i.e., H{x,y) — ${x) + *!/(?/), and as a 
byproduct we obtain a result for the case H{x,y) = A{x) + as well. We suppose in addition that 

Ii{h)=l fi(x)g(y)dx, for i = 0, 1, . . . , n - 1, 

where /o, /i, . . . , fn-i and g are analytic functions. (Note that the function depending on y is the same for all the 
1-forms. In the problems studied in the literature, the original family of Abelian integrals can be usually reduced 
to a family as above.) We will show that, in this case, some Chebyshev properties on fi and g (to be specified later 
on) transfer to li after the integration over the ovals. To fix notation, H is an analytic function in some open subset 
of the plane that has a local minimum at the origin. Then there exists a punctured neighbourhood V of the origin 
foliated by ovals 7^ C {H{x,y) — h}. We fix that H{0,0) — and then the set of ovals inside this, let us say, 
period annulus, can be parameterized by the energy levels h £ (0,/io) for some ho G (0, +00]. In what follows, we 
shall denote the projection of V on the s-axis by {xi,Xr). Similarly, {ye,yr) is the projection of P on the j/-axis. 

Theorem 1X1 is our main result and it applies in case that H{x,y) = $(a;) + '^{y). It is easy to verify that, under 
the above assumptions, x^'{x) > for any x £ (xe, Xr) \ {0} and y'i/'{y) > for any y € {ye,yr) \ {0}. Then ^ and 'i/ 
must have even multiplicity at 0. Thus, there exist two analytic involutions ai and (T2 such that 

$(a;) = <I>(cri(2:)) for all x e {xi,Xr) 

and 

^{y) = V['((J2(y)) for all y G {yt,yr). 

Recall that a mapping a is an involution it a o a = Id and a 7^ Id. Note that an involution is a diffeomorphism 
with a unique fixed point. In our situation we have that (7^(0) — 0. In what follows, given a function n, we define its 
balance with respect to a as 

iMcr{n){x) = K.{x) — k((7(x)). 

For example, if ct = —Id, then the balance of a function is twice its odd part. 

In the statement of Theorem m is related with the multiplicity of ^' at 1/ = 0. More concretely, we suppose 
that ^^(y) = ey^™ +0(1/^™) with e > 0. In addition, ECT-system stands for extended complete Chebyshev system in 
the sense of Mardesic [22], see Definition 12.11 for details. 

Theorem A. Let us consider the Abelian integrals 

Ii{h) = / fi{x)g(y)dx, i = 0,1, . . . ,n - I, 

where, for each h G (0, ho), 7^ is the oval surrounding the origin inside the level curve {<^{x) + ^(y) — h}. Let ai and 

(72 be the involutions associated to $ and VP, respectively. Setting go = g, we define gt+i = Then (Jo, Ji, . . . , In-i) 
IS an ECT-system on (0, ho) if the following hypothesis are satisfied: 

(a) (.^<Ti(|f),.^<Ti(|V),...,.?^ai(%i)) IS a CT-system on {0,Xr), and 

(6) (^^<,2(ffo),^<T2(ffi),---,^<T2(3n-i)) IS a CT-system on (0, and S3„J^go){y) = o(j/2'"("-2)). 
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To prove the result it is necessary to compute the derivative of each Abehan integral until order n — 1. The 
condition on ^0-2(50) (?/) at 1/ = ensures that the integral expression of this derivative is convergent, although it 
may be improper (see Remark 13.2 1) . Let us also point out that, since 0-2(2/) = + this condition is equivalent 
to require that g{y) - g{-y) = o(j/^™^""^>). 

Our second result deals with those Abelian integrals such that 

H{x,y) = A{x) + B{x)y'^"' and g{y) = y^"-'^ with s G N. 

Since H has a local minimum at the origin by assumption, B{0) > and A has a local minimum at 2; = 0. Thus, as 
before, there exists an involution a satisfying A{x) — A(^a{x)^ for all x £ {xi,Xr). 

Theorem B. Let us consider the Abelian integrals 

= fi{x)y^''^dx, i = 0,l,...,n-l, 

where, for each h £ (0, /lo), 7h is the oval surrounding the origin inside the level curve {A{x) + B{x)y'^"' = h}. Let g 
he the involution associated to A and we define 




Then (7o, /i, . . . , /n-i) is an ECT-system on (0, ho) if s > m{n — 2) and (fo, ^1, • ■ • , ^n-i) is a CT-system on {0,Xr). 

It is worth noting that although the condition s > m(n — 2) is not fulfilled in some situations, it is possible to 
obtain a new Abelian integral for which the corresponding s is large enough to verify the inequality. The procedure 
to obtain this new Abelian integral follows from the application of Lemma l4.ll We refer the reader to Example 14.21 
in which we explain in detail how to apply Lemma l4.1l to get a new Abelian integral with s > m(n — 2). 

The applicability of our criteria comes from the fact that the hypothesis requiring some functions to be a CT- 
system can be verified by computing Wronskians (see Lemma [2.3p . This simplifies a lot the problem of showing that 
a given collection of Abelian integrals has the Chebyshev property and in some cases it enables to reformulate the 
problem in a purely algebraic way (cf. Section |4|. 

In the literature there are a lot of papers dealing with zeros of Abelian integrals (see for instance O |6] |9l 1101 
1141 1231 124] and references there in) . In many cases, it is essential to show that a collection of Abelian integral has 
some kind of Chebyshev property. The techniques and arguments to tackle these problems are usually very long and 
highly non-trivial. For instance, in some papers (e.g. [4l[71[2l]) the authors study the geometrical properties of the 
so-called centroid curve using that it verifies a Riccati equation (which is itself deduced from a Picard-Fuchs system). 
In other papers (e.g. [51 1121 [TH] ). the authors use complex analysis and algebraic topology (analytic continuation, 
argument principle, monodromy, Picard-Lefschetz formula, ...). Certainly, the criterion that we present here can 
not be applied to all the situations (since the Abelian integrals need to have a specific structure) and, even in case 
that it is possible to apply it, sometimes the sufficient condition that we provide is not verified. However we want 
to stress that, when it works, it enables to extremely simplify the solution. To illustrate this fact, in Section |4] we 
reprove with our criterion the main results of three different papers. We are also convinced that this criterion will 
be useful to obtain new results on the issue. In this direction we tackle the program posed by Gautier, Gavrilov and 
Iliev [8] and we prove their conjecture in four new cases (see Subsection 14. ip . 

In several papers dealing with zeros of Abelian integrals (see [2] [3] |4l [21] for instance) , it is applied a criterion 
of Li and Zhang [19]. This criterion provides a sufficient condition for the monotonicity of the ratio of two Abelian 
integrals. In page 360 of the book of Arnold's problems [T], the criterion given in |19) is quoted as a useful tool that 
"despite its seemingly artificial form, it proves to be working in many independently arising particular cases". The 
translation of the result in [19] to the language of Chebyshev systems and Wronskians shows that it corresponds 
precisely to the case n = 2 of our criteria. Accordingly, using our formulation, their result becomes very natural: it 
shows that the Chebyshev properties of the functions in the 1-form are preserved after integration. In addition, as a 
generalization of their result, we hope that our criteria will be useful in many cases as well. Finally we remark that, 
although we suppose that the functions that we deal with are analytic, our results hold true for smooth functions 
with minor changes. 

The paper is organized as follows. Section [5] is devoted to introduce the definitions and the notation that we 
shall use. In particular we define the different types of Chebyshev property that we shall deal with and we establish 
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their equivalences with the continuous and discrete Wronskians (see Lemma [2.3|l . Theorems El and |B] are proved in 
Section El The main ingredient in the proof of Theorem |X1 is Proposition 13.31 that provides an integral expression 
for the Wronskian of a collection of Abelian integrals. Theorem [Bl follows as a corollary of Theorem |X1 Sectional is 
devoted to illustrate the application of our criteria. To this end, in Examples 14.21 IT3] and \4A\ we reprove the results 
of Iliev and Perko [8], Zhao, Liang and Lu 24] and Peng [21], respectively. Apart from showing the simplicity in 
the application of the criteria, our aim with these examples is twofold. First, to show that it is not necessary to 
know explicitly the involutions that appear in the statements. Second, to show that it is possible to reformulate the 
problem in such a way it suffices to check that some polynomials do not vanish. In Section |4l we also present some 
new results concerning the program of Gautier, Gavrilov and Iliev [^. Finally in the Appendix we give some details 
about the tools that are used in Section |4l namely, the notion of resultant between two polynomials and Sturm's 
Theorem. 



2 Chebyshev systems 



Definition 2.1 Let /o, /i, . . . , fn-i be analytic functions on an open interval L of R. 

(a) (/o, /i, . . . , fn-i) is a Chebyshev system (in short, T-system) on L if any nontrivial linear combination 

aofo{x) + aifi{x) + . . . + Q„_i/„_i(a;) 

has at most n — 1 isolated zeros on L. 

(b) (/o, /i, . . . , /n-i) is a complete Chebyshev system (in short, CT-system) on L if (/o, /i, . . . , fk-i) is a T-system 
for all fc = 1, 2, . . . , n. 

(c) (/o, /i, . . . , fn-i) is an extended complete Chebyshev system (in short, ECT-system) on L if, for all fc = 1, 2, . . . , n, 
any nontrivial linear combination 

aofo{x) + aifi{x) + . . . + au-ifk-i{x) 
has at most fc — 1 isolated zeros on L counted with multiplicities. 

(Let us mention that, in these abbreviations, "T" stands for Tchebycheff, which in some sources is the transcription 
of the Russian name Chebyshev.) □ 

It is clear that if (/o, /i, . . . , fn-i) is an ECT-system on L, then (/o, /i, . . . , fn-i) is a CT-system on L. However, 
the reverse implication is not true. 



Definition 2.2 Let /o, /i, . . . , fk-i be analytic functions on an open interval L of R. The continuous Wronskian of 
(/o,/i,---,/fc-i) at 2: G L is 



W[fo,fi,--- ,h-i]{x)=detlf^'\x) 



fo{x) 
fo{x) 



fk-i{x) 

f'k-l{x) 



The discrete Wronskian of (/o, /i, . . . , fk-i) at (2:0, a^i, . • . , Xk-i) G is 



D[,fo,.fi,--- ,fk-i]ixo,xi,...,Xk-i) = det(/j(a;i))(,^^ = 



/o(a;o) 
/o(a;i) 



fo{xk-i) 



h-i{xo) 
fk-iixi) 



fk-l{Xk-l) 



□ 



4 



For the sake of shortness, given any "letter" x and fc £ N we use the notation 

Xo,Xi, . . . ,Xk-l = Xk. 

Accordingly, we write 

W[fo,fi,--- ,h-i]{x) = W[f^]{x) 

and 

D[fo,fi,--- ,fk~i]{xo,xi,...,Xk-i) = -D[fk](xk) 
for the continuous and discrete Wronskian, respectively. The following result is well known (see [18ll22j for instance). 
Lemma 2.3. The following equivalences hold: 

(a) (/o, /i, . . . , /n-i) is a CT-system on L if, and only if, for each fc = 1, 2, . . . , n, 

D [fk] (xk) 7^ for all Xk G L'' such that Xi 7^ Xj for i 7^ j. 

(b) (/o, /i, . . . , /n-i) is an ECT-system on L if, and only if, for each = 1, 2, . . . , n, 

W [fk] {x) 7^ for all x£L. 

3 Proof of the main results 

The first part of this section is devoted to prove Theorem El Thus, unless we explicitly say the contrary, we suppose 
that H{x,y) = ^{x) + where ^'(y) = ej/^™ +0(1/^'") with e > 0, as mentioned before. Then, there exists a 

diffeomorphism /3 on {ye,yr) such that 

We take this diffeomorphism into account and we can write the involution associated to 'i' as 

a2{y)^rH-Piy))- 

In what follows, for each h € (0, /lo), we denote the projection of the oval -yn on the x-axis by {x'^,x'^). Therefore, 
Xi < x'^ < < x'l^ < Xr and $(a;^) = h. Moreover (see Figure [T]), if {x,y) G 7^, then 

y = Vhix) for ?/ > and y = y^ix) for y < 0, 

where 



yHx) := (± ''^2m{h-^x))'^ 



We note that y^^^x) = j/^ ((Ti(a;)) , where we recall that cti is the involution associated to $. We begin by the proof 
of the following result. 

Lemma 3.1. Let f and g be analytic functions on {xi,Xr) and [yi,yr), respectively, and let us consider 

I{h) = / f(x)g(y)dx. 



We set t{x):= f{x) — f(ai{xy)a[{x) and (^k~ -SSa^igk), where gk is recursively defined by means of gk+i = ^ with 
go=g. Then, */Co(y) = o(y2™("-2)), 

/('='(/.)= r\{x)ik{yt{x))dx /orfc = 0,l,...,n-l. 
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Figure 1: Notation related to the oval 7^. 



Proof. We prove the result by induction on k. We take the parameterization of the oval 7^ given by the mappings 
X I — > [x,yf^{x)), with the clockwise orientation, and we use y^{x) — 0-2 (3/^(2:)), to get that 

Hh) = fl' f{x)g{y^{x))dx+ f{x)g{y+{x))dx^ f{x) {g{y) - g{a2{y)))\^^^+^^^dx 
= J_f{x) {g{y) - g{'^2{v)))\y^y+i^^^dx + f{x) {g{y) - g{o2{y)))\y^^+^^^dx 

^f{'^^(u))(j[[u) {g{y) - g{(T2{y)))\y^y+i^^^i^^^^du + f(x) {g{y) - g{cr2{y)))\y^y+^^^dx, 

where in the last equality we performed the change of variable x = ai{u). Thus, since 2/^(o"i(m)) = y^{u), the above 
expression yields to 

Hh) = 1^^' im - f{ar{x))a[{x)) {g{y) - g{a2{y)))\^^^+^^^ dx = J^^" £{x).<M,, (<?) {yUx))dx. 

This expression proves the result for k — 0. We assume now that the result holds true for k < n — 1. On account of 
the hypothesis about the order of ^0 at y = 0, an easy computation shows that ^k{y) = >^(T2 (fffc) (y) = o(j/^™'^"~'^~'°-') . 
The fact that 2m(n — 2 — fc) ^ enables us to differentiate the expression of and we obtain 



d 



dx 



(Let us note that in the second equality we use that y^{x) = at x = a;^ because ^{x'^) — h and 'I'(i/^(a;)) = h for 
all h.) Finally, since 

^Uy) = gUy) - 9'ki<^2{y))a'2{y) = a'^iv) ^ M'^^iv))^^^^ = = 

the result for fc + 1 follows and the proof is completed. ■ 
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Remark 3.2 It is worth making some comments on the expression of the (n — 1) derivative of I{h) given by 
Lemma [3.11 The condition 3§cT2{9o){y) = (,o{y) ~ o(j/^'"'"~'^-') guarantees that the integral 



r(n-l) 



[h) 



despite it may be improper, is convergent. Indeed, by this condition, the Taylor series of at y = begins at least 
with order 2m(n - 2) + 1, i.e. ^o{y) = Ai/^'"'""^^+^ + . . . with A / 0. To construct gk+i{y), we derive gk{y) and 
divide it by ^''(y), which vanishes at i/ = with multiplicity 2m — 1. Hence, it turns out that ^n-i = i^o-2(5n-i) is 
not analytic at y = but meromorphic. However, due to the mentioned condition, the pole has at most order 2m — 1. 
We note that y^{x) = at x = x^ because ^{x'^) — h. More precisely, we take 7^ also into account and it 

is easy to show that 



1™, 2™r 

X — y ^ 



Accordingly, although ^n-i (j/^ (*)) i^^y tend to infinity as x — > x^, the derivative ^Hh) is given by a convergent 
integral. □ 



Let us consider now 



Ikih) = / fk{x)g{y)dx, for fc = 0, 1, . . . , n 



where g is an analytic function on {ye,yr) and each fk is an analytic function on (xi,Xr). The next result provides 
an expression of the Wronskian of (/o, /i, . . . , In its statement, is defined as in Lemma l3.ll i.e. we set 

fli+i = ^ with go = g, and := ^^sCffi)- Moreover 

Afe(/i):= {xk G K''^ < a;o < < . . . < Xk^i < x^}. 

Proposition 3.3. Let us assume that ■dSa2{g){y) = o(j/^'"'-"~^-') . Then, for each fc = 1,2, ...,n, the Wronskian of 
(/(), /i, • . • , h-i) at h £ (0, h(i) is given by 



W[l^]{h) 



D [fk] (xk)I> kk] (yk) dxo dxi--- dxk-i, 



Afc{h) 



where yi = y'^{xi) and £i{x) = fi{x) - ft[ai{x))(T[(x). 

Proof. Fix k £ {1, 2, . . . , n} and let Sk be the symmetric group of k elements. We take the definition of determinant 
into account and we apply Lemma l3.1l to show that 



W[l^]{h)^Aet[lf\h)) = E«gnWn4(l)W 

= E '^g'^i(^)n / ^T{i){x)£,i{yl{x))dx 



=0 
fc-1 



= E sgn(T)|| / iT(i:){xi)(,^{yl{x^))dxi 



^i{yi) dxo dxi 



■ dxk 



-D[fk] (xk) Y\ iiiVi) dxo dxi ■ ■ ■ dxk 



'I0,x + ]- 

At this point, for each permutation t £ Sfc we define ipr '■ K'^ — > K'^ as 

TprixoyXi, . . . ,Xk-i) = (a;^(o), a;^(i), • • • ,Xr(k-i)), 
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which is clearly an invertible mapping. We note that 

[o,x+]'=\7^= y v.(Afe(/.)), 

where 7?. is a subset of R'' with Lebesgue measure equal to zero. Accordingly 



W[l^]{h)= f ■■■ [ D[l^]{^^)W(,,{y,)dxodx^---dxk-i 

= /■■■/ D[ei^](xi^)T^(,i{yi)dxodxi ■ ■ ■ dxk-i- 

Next, in each integral of the above summation we perform the coordinate transformation Xk = ^/jrCuk) (i.e., Xi = Ur(i) 
for i = 0, 1, . . . , fc — 1), so that 

where Vi = Vhi^i)- (Here we use that the absolute value of the determinant of the Jacobian of tpr is identically 
one.) Finally, we remark that Of^k] (i/'T(uk)) ~ sgn(T)_D [fk] (uk) and we take the properties of the determinant 
into account to prove that 

W [ik] (^) = ^ I " I sgn(r)_D [^k] (uk) W £,i {vr^i)) duo dui ■ ■ ■ duk-i 

= /••■/ D [^k] (uk) y2 sgn(r) TT («T(i)) '^"o dm ■ ■ ■ duk-i 

= / ■•■ / D[^k] (uk)-D[Ck](vk) dMorful • ■ • dufc_i, 
and this last identity proves the result. ■ 

Proof of Theorem |A1 We claim that the assumptions (a) and (b) imply that the Wronskians W [ik] (h) for k = 
1,2, ...,n are different from zero at any h £ (0, fto)- On account of (b) in Lemma 12.31 this fact will prove that 
(Jo, /i, . . . , In-i) is an ECT-system on (0, ho). 
From Proposition 13.31 

W[l^]{h)^ [■ ■[ D[ei,]{x^)D[myi,)dxodx,---dxk-i, 

where recaU that yi — y'^ixi) — f3^^ ^'^2m{h — <^{xi))^ . On the other hand, x i — > f5~^ ^'^2m(/i — (^{xi))^ is 
decreasing on (0, Xr) and, therefore, in the above integral we have that 

< xo < a;i < . . . < xu-i < and < yk-\ < yk-2 < ■ ■ ■ < yo < Vh ■ 

We note at this point that £i(x) = ^' {x)£i§a-i{-^) (x) because 

£,ix) = ,n{x) - ,n{a,ix))a{{x) = Mx) - f,{a^ix)) = $'(x) ((A)(x) - (A)(ai(x))). 

Since $'(a;) 7^ for any x G {xe,Xr) and, by assumption, (^^o-i(^) , .^<ti(^) , • • . , >^cti(^%7-^)^ is a CT-system on 

{0,Xr), so it is {£o,£i, ■ ■ ■ ,£n-i)- The second assumption ensures that (^0, ^1, . . . , ^n-i) is a CT-system on (0,2/, ) 
because, by definition, = ^0-2 (gi). Therefore, we apply statement (a) in Lemma 12.31 and it turns out that 

D[£^] (Xk)£i[a] (yk) / for aU Xk G Ak{h). 

Since Ak{h) is connected, we have shown that [ik] (h) 7^ and the result follows. ■ 
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Proof of Theorem IbI This result is in fact a corollary of Theorem 1X1 We note that B{x) > for x £ {xe,Xr). 
Thus the coordinate transformation (m, v) = x(^j v) ■= (^i ^'\/2mi3(a;) y) is well defined and verifies eu := X~^(7h) 
{A(u) + j^u^™ = h). Accordingly 

Hh)^ j Hx)y^"-'dx = {2m)^ I (^t^\u)v^'-^du. 

Following the obvious notation, we can apply Theorem 1^ with 

% = ^i=T, g[v)=v^''\ ^ = A, <l>{v) = ^v^^, ai = a and <72 = -7d. 

B 2m 

Clearly the hypothesis (a) in Theorem|X]is guaranteed by the assumption on (.i = 38,^ ('^^ • Let us turn now to the 

hypothesis (fe). We take (72 = —Id and — u^'""^ into account and one can easily show that 'gi{v) = Ci-i;^'-''"""'^^ 

for some positive constant Ci, so that S§a2{{ji){'") = 2ci77^'°~""^~"^ . Hence, ^^§'^2(50), •^^CT2(fli)5 • • • , •^CT2(5n-i)^ is 

clearly a CT-system on (0,+oo). Since the condition s > m{n — 2) implies that ^ct2(?)(^) = 2^^""^ = o(u^'"*"~^'), 
the hypothesis (6) in Theorem |X] is satisfied as well. Therefore, we apply Theorem |X] and we can assert that 
(Jo, /i, . . . , In-i) is an ECT-system on (0, ho) as desired. ■ 



4 Applications 

The following lemma establishes a formula to write the integrand of an Abelian integral so as to be suitable to apply 
our results. 

Lemma 4.1. Let jn be an oval inside the level curve {A{x) + B{x)y^ = h} and we consider a function F such that 
F/A' is analytic at x = 0. Then, for any fc G N, 

/ F{x)y''~^dx = / G{x)y''dx 
J ih J ih 

where G(x) = l{^)'{x)-{S^){x). 

Proof. If {x,y) G 7h C {A{x) + B{x)y^ = h} then £ = -^^±^^S2^^ and accordingly 

d{g{x)y^) =g'{x)y''dx + kg[x)y^~^dy 

^(ff'W - y''dx-^{^){x)y''-'dx. 

We take F{x) = ^{^){x) in the above equality, we use that J d(cj{x)y''^ = and the result follows. ■ 

From now on we shall often compute the resultant between two polynomials and we shall apply Sturm's Theorem 
to study the number of roots of a polynomial in an interval. The interested reader is referred to the Appendix for 
details. 



Example 4.2 Iliev and Perko study in llj symmetric Hamiltonian systems perturbed asymmetrically. More con- 
cretely, systems of the form 

x = y, 

y = ±{x ± x-'') + \iy + \2X^ + \j.xy + XiX^y, 

where \j{e) = 0(e), and they prove that at most two limit cycles bifurcate for small e 7^ from any period annulus 
of the unperturbed system. There are three different cases to consider depending on the phase portrait of the 
unperturbed system: the global center, the truncated pendulum and the Duffing oscillator. This latter case gives 
rise to two different types of period annuli (see Figure [^jl . In this example we study the so-called interior Duffing 
oscillator. Theorem 1.3 in shows that at most two limit cycles bifurcate from either one of the interior period 
annuli. 
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Figure 2: The period annuli in the Duffing oscillator. 



If we perform a translation to bring the center on the right half-plane to the origin, the Hamiltonian function of 
the unperturbed system becomes 

H{x,y) = A{x) + B{x)y^ with A{x) = x'^ + \x'^ + x^ and B{x) = i. 

The projection of the period annulus of this center is (— f , \/2 — f) and ha — A{—1) = f /4. 

From Theorem 2.1 in [11], it follows that the first non-identically zero Melnikov function is a linear combination 
of Ii{h) = J^^x^ydx for i = 0,1,2. Thus, Theorem 1.3 in [TJ will follow if we prove that {/o,/i,/2} is an ECT- 
system. Additionally, this fact implies that there are values of the parameters for which exactly 0, 1 or 2 limit cycles 
bifurcate from the period annulus. To this end we will apply TheoremlB] but we note that in this case m = 1, n = 3 
and s = 1, so that the hypothesis s > m(n — 2) is not satisfied. This is easy to overcome because 

Io{h)= f ydx^\ [ {A{x) + B{x)y'^)ydx= \ [ A{x)ydx+\[ B{x)y^dx, 

and then, we apply Lemma |4 . 1 1 with fc = 3 and F = ^4 to the first integral above, to get 

~ If x^ + 2x + 2 3 , ^ 1 /■ 1 3 , 1 /■ f ^ ^ 3 , f I ^ 7x' + Ux + ^ 

Io(h) = T -TTTr — r-^y "^+r / T^y dx= - fo{x)y dx with fo{x): 



(It is not possible to apply Lemma [4.11 directlv to Iq because then we must take F = 1, and in this case F/A' is not 
analytic at a: = 0.) Exactly in the same way we obtain 

^hj -^'^"^^^ ./i(^):= I2{x + iy ' 

7 1 /■ . . N 3 , . ^ ^ a:'(to'+20a; + 12) 
h(h) = T / f2(x)y dx with ,f2(x):= — ^ - 



hj.J^''' 12(x + l)2 • 

We set Ii{h) = J fi{x)y^dx and it is clear that {Iq, Ii, I2} is an ECT-system on (0,/io) if and only if so it is 
{Iq, Ii, 12}- We can now apply Theorem iBl because s = 2 and the condition s > m{n— 2) holds. Thus, setting 



we have to check that {£0,^1, £2} is a CT-system on (O, y/2 — l). Here a is the involution associated to A and we 
used that B is constant. (In this example we can compute the involution explicitly but we do not use it because we 
want to show that it is not necessary to apply our result.) As a matter of fact we will show that {^0,^1,^2} is an 
ECT-system because a continuous Wronskian is easy to study. In order to compute the three Wronskians, we write 
£i{x) = Li(^x,a{x)) with Li{x,z) = {^)ix) — (^^){z). Moreover, due to 

A(x) - A{z) = l{x - z){x + 2 + z){x^ + 2x + 2z + z^), 



10 



it turns out that z — a{x) is defined by means of q{x, z) := x'^ + 2x + 2z + z'^ — 0. Accordingly, since o-'{x) — — f^, 
we have that W[£i] (x) = LVi (a;, o"(x)) with uii {x, z) being a rational function for i — 1,2,3. The resultant with respect 
to z between q{x, z) and the numerator of uj^ix, z) is rz{x) — &A:X^'^{x + 2Y^p'i{x) with 

p-i{x) = 441 a;^° +8820 + 79380 a;^® + 423360 x^^ + 1481685 x^*^ + 3555024 a;^^ + 5918640 a;" 
+ 6740160 + 4976155 x^'^ + 1881540 - 892716 - 3303200 x^ - 4779945 
- 3240840 x^ + 601960 x" + 2523360 x^ + 1158080 x" - 414400 x^ - 414400 x^ + 44800, 

and by applying Sturm's Theorem we can assert that p-i{x) 7^ for all x G (O, \/2 — l). Thus, ll1z{x,z) = and 
q{x,z) = have no common roots, and this fact implies that W[^3](x) 7^ for all x € [0,y/2— l). The resultant 
with respect to z between q{x, z) and the numerator of uj2{x, z) is r2(x) = 32x^{x + 2)'^j32(x) with 

P2(x) = 49 x^^ + 588 x" + 2940 x^" + 7840 x^ + 11650 x** + 8528 x'^ 
+ 496 x^ - 3520 x^ - 1915 x'* - 620 x^ - 620 x^ + 360, 

and using Sturm's Theorem it follows that p2 does not vanish on (0, \/2 — 1). Exactly as before, this fact shows that 
VF[i'2](x) 7^ for all X £ (O, \/2 — l)- Finally, the resultant with respect to z between q{x,z) and the numerator of 
a;i(x,z) is 

n (x) = 2 x^ (x + 2)^ (49 X* + 392 x'^ + 1176 x*^ + 1568 x^ + 659 x* - 500 x^ - 500 x^ + 80) 

and, thanks to Sturm's Theorem again, we can assert that it does not vanish on (0, \/2 — 1). This proves that 
W[^i](x) = £o(x) / for aU x G (0,^- 1). Consequently {4,^1,^2} is an ECT-system on (0,^2- 1) and by 
applying Theorem [B] {Io,Ii,l2} is an ECT-system on (0, 1/4). Therefore, the first Melnikov function has at most 
two zeros counting multiplicities. □ 

Example 4.3 Zhao, Liang and Lu study in [24] the system of planar differential equations 

x = 2xy + £(^Y^ ay(e)x'y^^, 

2) = 6x - 6x^ - + ^ bij{e)x'y^y 

i+j^2 



The unperturbed system (i.e., with e = 0) has a center at (1,0) whose period annulus is bounded by a cuspidal loop 
and they prove (see Theorem 1.2 in [24 j ) that the maximum number of limit cycles emerging from its period annulus 
for e ~ is two. 

Our goal is to reobtain this result by applying Theorem |B] To this end, we bring the center to the origin by 
means of a translation, so that the unperturbed system is Hamiltonian with 

H{x,y) = A{x) + B{x)y'^, where A{x) = x^(3 + 2x) and B{x) = x + 1. 

The projection of the period annulus is now (—1, 1/2) and the energy level of the polycycle in its outer boundary is 
ho = A{—1) = 1. By Theorem 3 in 16 , the upper bound for the number of limit cycles is equal to the maximum 
number of zeros for h £ (0, 1), counted with multiplicities, of any non-trivial linear combination of 



Ii{h) = / (x + if'^ydx for i = 0, 1, 2. 
•'Jh 

Accordingly, the result in [24] will follow once we show that {/o,/i,/2} is an ECT-system on (0, 1). By applying 
Lemma [4. II the same straightforward manipulation as before shows that Ii{h) = j^Ii{h) where 

U{h) = f Mx)y'dx 

with 

16x2+35x-H24 20x2 41x 24 , 2 , , 

^"^""^^ {x + lf ' ^TTl and /2(x) = 24X'' +47x-h24. 
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It is clear that {lo, Ii, I2} is an ECT-system on the interval (0, 1) if, and only if, so it is {7o,/i,/2}- On account of 
Theorem IB] this will follow once we check that {^0,^1,^2} is an ECT-system on (0, 1/2), where £i = SS^ ^ a'b^/^ ) ' 
Note that A{x) — A{z) = (x — z)(2x'^ + 2zx + 3x + 2z^ + 3^), so that z = cj{x) is implicitly defined by means of 
q{x, z) ■- 2x^ + 2zx + ix + 2z^ + 3z = 0. Thus 

,, , dz 4x + 2z + 3 
a (x) = — = . 

^ ^ dx 'iz + 2x + 3 

Taking this into account, some computations show that, for i — 1, 2, 3, VK[i'i](2:) = LUi(^x, o-(x)) with uji{x, z) being a 
rational function of ii = y/x + 1 and v — \/z + 1, say Ri(u,v). Note that x 1 — > y/x + 1 maps (0, 1/2) to (1, ^3/2). 
The resultant with respect to v between the numerator of Ri{u, v) and q{u'^ — l,v^ — 1) is a polynomial ri{u) that, by 
applying Sturm's Theorem, has no roots on (1, ^3/2). (For the sake of shortness we do not give here the expression 
of these polynomials.) Hence, it is proved that does not vanish on (0, 1/2) for i = 1, 2, 3. By Theorem IB] this 

reasoning proves the mentioned result of Zhao, Liang and Lu. □ 

Example 4.4 Peng studies in i21.. the system of planar differential equations 

r = — y — Sx'^ — + £{fJ-ix + fJ.2Xy), 



y = x{l- 2y) + £^.zx^. 

The unperturbed system (i.e. when e = 0) has a center at the origin and the author proves (see Theorem A in [21]) 
that two is the maximal number of limit cycles which bifurcate from its period annulus for e ~ and that there 
are perturbations with exactly 0, 1 or 2 limit cycles. To this end, he first shows that by means of the projective 
coordinate transformation (x,y) t-^ (74:2' 2[x+2) ) ^ non-constant rescaling of time the above system reads for 

x = 2{x + 2)y + e^.z{x + 2)y^ 

y = -x - fx^ - 2/2 _^ e (^ii {x + 2) + fx + ^liy^). 

The unperturbed system is now Hamiltonian with a center at the origin whose period annulus is bounded by a saddle 
loop. We have written the transformations so as to directly apply Theorem [B] The Hamiltonian function of the 
unperturbed system is 

H{x,y) = A{x) + B{x)y'^ with A{x) = \x^{x + 2) and B{x) = x + 2. 

The projection of the period annulus is (—4/3,2/3) and the polycycle at its outer boundary has energy level ho = 
A(2/3) — 8/27. It is very easy to show that the first Melnikov function is a linear combination of 

Ii{h)= {x + 2yydx ior i = Q,l,2. 



Hence, the aforementioned result will follow once we check that {Jo,/i,/2} is an ECT-system on (0,/io). By using 
Lemma [4.11 exactly as before, Uiji) — j^Ii{h) where 

Hh) = / Mx)y''dx 
•'ih 

■4--U J- / \ 2{a: + 2){15a:^-|-42a;-|-32) r / n 4(3: + 2)^ (9a:^ +23a;-|-16) j r / \ 2(a;-|-2)^ (21a:^ -|-50a;-|-32) ^ ■ l7 7 7 ^ 

with/o(a;) = ^ 3(3:^+4)^ ' /n^) = 3(3:r+4)a and/2(x) = ^ 3(3^+4)^ ■ Onceagam, {7o,/i,/2j 

is an ECT-system on (0, ho) if, and only if, so it is {/o, Ii, h}- The involution associated to A is z = o-{x) given by 
q{x, z) := x'^ + xz + 2x + z'^ + 2z — because A{x) — A{z) — j{x — z)q{x, z). Thus 

, dz z + 2x + 2 
^ ' dx x + 2z + 2 

and, setting it — SSa (^ ^,^3/2 ^ i we have to verify that W[£i] does not vanish on (0, 2/3) for i = 1, 2, 3. It can be shown 
that, for i = 1,2,3, M^[i?i](a;) — u)i(^x,a{x)) with LOi{x,z) being a rational function of u = y^x + 2 and v = \^z + 2, 
say Ri{u, v). We note that x 1 — > \Jx -\-2 maps (0, 2/3) to {\p2, s/%J?>). The resultant with respect to v between the 
numerator of Ri{u, v) and q{u^ — 2,v'^ — 2) is a polynomial ri{u) that, by applying Sturm's Theorem, has no roots on 
(\/2, \/8/3). Therefore, W[fi] does not vanish on (0,2/3) for i = 1,2,3. By Theorem [B1 we have proved the result 
of Peng in pT]. □ 
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4.1 Results on the program of Gautier, Gavrilov and Iliev 



Our last examples of application come from the paper of Gautier, Gavrilov and Iliev where a program for finding 
the cyclicity of the period annuli of quadratic systems with centers of genus one is presented. They give a list of 
the essential perturbations of these centers (i.e., the one-parameter perturbations that produce the maximal number 
of limit cycles), together with the corresponding generating function of limit cycles (i.e., the Poincare-Pontryagin- 
Melnikov function). Since some cases have been already solved in the literature about the problem, this list includes 
only the open cases, a total of 26. They conjecture that the cyclicity of these period annuli is two, except for some 
particular cases in which it is three (cf. Conjecture 1 in page 12 and Conjecture 2 in page 17). In their Theorem 3, 
two quadratic reversible systems with a center are considered, denoted by (rll) and (rl8) in the list, and they show 
that, in both cases, the upper bound of the number of limit cycles produced by the period annulus under quadratic 
perturbations is equal to two. We are going to reobtain this result for the case (rll) by using our criterion. Moreover, 
we prove their conjecture in four new cases in their list, namely (r7-rl4), (rl5), (rl7) and (rlv3). In fact. Theorem 
[B1 is likely to be applied in many of their cases but we have only been able to directly show that the functions on 
the integrand satisfy the Chebyshev condition in the five mentioned cases. We remark that our criterion gives a 
sufficient condition for the Abelian integrals to be an ECT-system. 

Case (rll) We translate the center to the origin, so that the first integral of the unperturbed system is 

H{x,y) = A{x)+B{x)y'' with A{x) = and B{x) = ^j^. 

They show that the cyclicity of the period annulus under quadratic perturbations is two. This will follow once we 
show that {/o,/i,/2} is an ECT-system on (0, 1/6), where 

h(h) = I [x + iy-'^ydx. 

The projection of the period annulus of the center at the origin is ( — 1/3, +00). By applying Lemma [4. II once again, 
Ii{h) = ^h{h) where Uih) = J^^ fi{x)y'^dx with 

^ , , 5x^ + 13x + 24 ^ , , 7a;^ + 19x 4- 24 , ^ , ^ 9x^ + 25x + 24 

It is clear then that it suffices to show that {/o,/i,/2} is an ECT-system on (0, 1/6). With this aim in view, let us 
note that A{x) — A{z) = Q(x+iyi(z+iy^ with q{x, z) := 3x^z + x'^ + lOxz -|- 3x -I- 3xz^ + z'^ + 3z, so that the involution 
z = a{x) associated to A satisfies q(^x,a{x)^ = 0. Taking this into account, we get that 

^ _ dz _ x{z + 

As before we must compute the Wronskians M^[£i](a;) for i = 1,2,3, where ii = SSa ^ j^i ^3/2 j 1 ^-nd then show that 

they do not vanish for x G (0,+C)o). In this case VF[i?i](a::) = u!i(x,a{xyj with LUi{x,z) being a rational function of 
u — x/aT+T and v — \/ z + 1, say Ri{u,v). The resultant with respect to v between the numerator of Ri{u,v) and 
q{u^ — l,v^ — 1) is a polynomial ri{u). Since the mapping x 1 — > ^^x + \ sends (0, -l-oo) to (l,-|-oo), the result will 
follow once we show that these polynomials ri(it) do not vanish on (l,+oo). This latter fact is deduced from the 
application of Sturm's Theorem. □ 

Let us mention that we have studied the case (rl8) as well (the other case that contemplates Theorem 3 in [8]), 
but it seems that it cannot be solved by using the criterion given by our Theorem [B] Of course, the success in the 
application of this criterion depends on the particular problem studied, but we want to stress that, when it works, it 
enables to extremely simplify the solution. For instance, the proof of Theorem 3 takes eight pages of highly nontrivial 
arguments. From now on, for the sake of brevity in the exposition, we omit many of the explanations on the way to 
apply our criterion since they are a verbatim repetition of the previous examples. 

Cases (r7-rl4) and (rl5) The first integral is shared by the two cases and, after we translate the center at the 
origin, it reads for 

H(.,,) = ^ + ^'(3.^ + 8x + 6) 



2 12 
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The cyclicity of the period annulus, whose projection on the a;-axis is the interval (—1, 1/3), is two if we prove that 
{Jo, 7i, I2} is an ECT-system for h G (0, 1/12), where 

Ii{h) = / {x + iy^^ydx for the case (r7-rl4), 

■'7(1 

Uih)^ f ix + iy-^ydx for the case (rl5). 

We apply Lemma |4. II to the Abelian integrals given by Ii(h) — hli(h) in order to write them in the form Ii{h) = 
f fi{^)y^dx. We have that: 



f rr, ^ , f 211^ + 63x2 + 64x + 24 

L = L — ^^(^' — ' ' 

/ Hix,y)y^=[ (2^ + 3)(9x2 + 14x + 8) 3 
A. ^ + ^ l„ 36(x + l)4 



, , 15x^ + 47^2 + 52a; + 24 3^ 

J^(x,y) — - — I — y ax, 

V 'W^ , 1^2 / 36(x + l)5 



/ 



ydx 


f (2: 


x + 1 




ydx 


-j 


(X+1)2 


■'ih 


ydx 




(x + iy^ 


•llh 


ydx 


-I 


(x + l)4 





, , 122:^ + 39^2 + 46a; + 24 3, 

Hi-,y)7rT^-l S6ix + ir 

H(x y) - I to^ + 31^^+40x + 24 3^ 

^(^^y) ,^ , 1^4 - / 36(a; + l)7 ^ 

Some computations show that the involution a defined by A{x) :— H{x,0) satisfies g(a;,(j(a;)) — with q(x,z) :— 
Zz^ + "ixz^ + Sz^ +?,x^z + 8a:z + &z + 3x^ + Sa;^ + 6x. We use resultants and Sturm's Theorem in order to check that 
the corresponding Wronskians have no zeros on the interval (0, 1/3). □ 

Case (rl7) Once the center is translated to the origin, the first integral reads for 

iJ(a;,y) = ^ + i^^±M. 
\ 2 6 

Setting Ii{h) — J^^ix + iy~'^ydx, the cyclicity of its period annulus is two if we prove that {lo, Ii, I2} is an ECT- 
system on (0, 1/6). By Lemma l4.ll we have that Ii{h) = J^^ fi{x)y'^dx, with 

, , , 5a;2 + 13a; + 12 ^ , , 7a;2 + 16a; + 12 , , , , 9x^ + 19a; + 12 
/oW = (, + 1)5 ' ^(-) - (a; + 1)4 ^'^d ^(-) - (. + 1)3 • 

In this case, the involution a defined by A{x) := H{x,0) satisfies g(a;,cr(a;)) — where q{x,z) ~ 2z^ + 2xz + 3z + 
2x^ + 3a;. The projection of the period annulus on the a;-axis is (—1, 1/2) and, thus, we are done if we show that the 

functions £i = £Scr ^ j^,g3/2 ^ form an ECT-system in (0, 1/2). Once again, the involution can be explicitly written, 
but we prefer to use resultants and Sturm's Theorem because it provides an algebraic procedure to check that the 
Wronskians do not vanish on (0, 1/2) for i = 1, 2, 3. The proof of this fact is omitted for the sake of shortness. □ 

Case (rlv3) After the center is translated to the origin, the first integral becomes 

H{x,y)=x\2-x^) + ^{l + xyy\ 

Since A{x) := H{x,0) is an even function, we have that (t(x) = —x and this simplifies a lot the computations. 
The projection of the period annulus on the x-axis is (—1, 1). In order to prove that its cyclicity under quadratic 
perturbations is two, we are lead to show that {/o,/i,/2} form an ECT-system for h G (0,1), where Ii{h) = 
/^^ Mx)y'dx with fo{x) = (5.--2.3-9.^+4.+8)(.+i) ^ ^ (r.^-i3._-+8)(.+i) ^ '^-^^^ (;l\-^2-+^ To 

this end, by applying Theorem IbI and taking <j{x) — —x into account, it suffices to show that the functions 

„ , , 5x^-8x^ + 7x^ + 8 „ , , 7x'' - ISx^ +8 , „ , s 5x'^ - 9x^ + 8 
^o(a;) = ——, — --ni"' ^i(^) = ~f TW? — r~iv7 ^^'^ ^2(^) = 



2x(x-i)5(x-hi)5 ' ' x{x-iyix + iy-' ^' x{x-iY{x + iY 
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form an ECT-system on (0, 1). It is easy to see that li does not vanish on (0, 1). The Wronskian associated to l\ 
and li is the rational function 

"^y^-M^^^ = i8.(x-l)B(x+l)B ' 

which has no zero on (0, 1) by virtue of Sturm's Theorem. Finally 

2(512 - 1632a;^ + 2200x* - 1617x^ + 693x* - 175a;^° + 352;^^) 



9(a;-l)iB(a; + l)iB 



which neither vanishes on (0,1), again by using Sturm's Theorem. As desired, this shows that (-£2,^1,^0) is an 
ECT-system on (0,1). □ 



5 Appendix 

5.1 Resultant of two polynomials 

Given two polynomials p, g G C[a;,2/], say 

p(x-) — aa'j-^ + arz^'^ + . . . + am, with ao / 0, 
(j(x-) = 6().x-" + fei.x""-' + . . . + 6„, with &o / 0, 

where a,, 6, £ C[j/], the resultant of p and q with respect to x, denoted by Res(p, g, x) is the (m + n) x (m + n) 
determinant 



Res(p, q, x) = det 



ai ao 
02 ai 
: a2 



60 

62 61 

ao : 

ai 6„ : 

a2 6n 



V 



60 

61 
62 



&n / 



where the blank spaces are filled with zeros. The three basic properties of the resultant are: 

1. Kes{p,q,x) is an integer polynomial in the coefficients of p and q. 

2. Res(p, q,x) = if, and only if, p and q have a nontrivial common factor in C[x, y]. 

3. There are polynomials A,B € C[x, y] such that Ap + Bq = Res(p, q, x). Moreover the coefficients of A and B 
are integer polynomials in the coefficients of p and q. 

Resultants can be used to eliminate variables from systems of polynomial equations. As an example, let us 
suppose that we want to study the following system of two polynomial equations with two variables: 



xy ■ 



x^ +y 



1 = 0, 

2 4: 



0. 



Here we have two variables to work with, but if we regard p(a;,j/) := xy—1 and q{x, y) := x^+y^—4: as polynomials in x 
whose coefficients are polynomials in y, we can compute the resultant with respect to x to obtain Res(p, q, x) = y* — 

4y^ + 1. By the third property above, there arc polynomials A,B £ C[x, y] such that A{x, y)p{x, y) + B{x, y)q{x, y) — 
y'^—Ay^ + 1. Accordingly, y'^—Ay^ + 1 vanishes at any common solution of p = g = 0. Thus, we can solve — + 1 = 
and find the y-coordinates of these solutions. 
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5.2 Sturm's Theorem 

A sequence {/o, /i, . . . , fm} of continuous real functions on [a, b] is called a Sturm's sequence for f — fo on [a, b] if 
the following is verified: 

1. /o is differentiable on [a,b]. 

2. fm does not vanish on [a,b]. 

3. If /(2;o) = with xo € [a, 6] then fi{xo),fo{xo) > 0. 

4. If /i(a;o) = with xq G [a, 6] then /i+i(a;o),f,_i(a;o) < 0. 

Sturm's Theorem. Let {/o, /i, • ■ • , /m} be a Sturm's sequence for f = fo on [a,b] with f{a)f{b) ^ 0. Then the 
number of roots of f on (a, b) is equal to V{a) — V{b), where V{c) is the number of changes of sign in the sequence 
{/o(c),/i(c),...,/™(c)}. 

There is a simple procedure to construct a Sturm's sequence in case that / is polynomial. Indeed, if p{x) is a 
polynomial of degree n, we define the sequence {po,Pi, • ■ • ,Pm} with m ^ n in the following way. We set po = p, 
Pi = p' and 

Pi-i{x) = qi(x)pi(x) -pi+x{x), for i = 1,2, . . . ,m - 1, 

Pm~l{x) = qm{,x)pm{,x), 

where qi{x) and pi+i{x) are the quotient and the remainder (the latter with the sign changed) of the division of 
Pi-i{x) by Pi{x), respectively. The construction of this sequence ends when the remainder is zero, i.e., Pm+i = 0. 
In this case, since this is essentially Euclides' algorithm, pm is the greatest common divisor of po and p\. If all the 
zeros of p are simple then pm does not vanish and it is easy to show that {po,Pi, • • • ,Pm} is a Sturm's sequence for 
p on any interval. If p has zeros with multiplicity then vanishes. Since pm divides po and pi, it also divides pi 
for i = 2, 3, . . . , m. In this case, we set pi = pi/pm and it follows that {po,pi, . . . ,Pm} is a Sturm's sequence for p on 
any interval. 
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